Line-wave waveguide engineering using Hermitian and non-Hermitian metasurfaces

Line waves (LWs) refer to confined edge modes that propagate along the interface of dual electromagnetic metasurfaces while maintaining mirror reflection symmetries. Previous research has both theoretically and experimentally investigated these waves, revealing their presence in the microwave and terahertz frequency ranges. In addition, a comprehensive exploration has been conducted on the implementation of non-Hermitian LWs by establishing the parity-time symmetry. This study introduces a cutting-edge dual-band line-wave waveguide, enabling the realization of LWs within the terahertz and infrared spectrums. Our work is centered around analyzing the functionalities of existing applications of LWs within a specific field. In addition, a novel non-Hermitian platform is proposed. We address feasible practical implementations of non-Hermitian LWs by placing a graphene-based metasurface on an epsilon-near-zero material. This study delves into the advantages of the proposed framework compared to previously examined structures, involving both analytical and numerical examinations of how these waves propagate and the underlying physical mechanisms.

advantage in terms of singular field enhancement.This characteristic, where infinite energy is concentrated along a line, represents another crucial aspect of these modes [21][22][23][24][25] .
In this context, we introduce a graphene patch design to present LWs in the terahertz and infrared ranges.While the implementation of LWs is typically limited to the terahertz range, we expand the operational range of these waves into the infrared domain by designing a dual-band line-wave waveguide.Since the classical metal conductivity model is no longer valid at higher frequencies, we use a graphene-based metasurface in an effort to implement the structure.Our designs consist of a graphene metasurface positioned atop two distinct substrates.One substrate is an ENZ substrate, which is a type of material with a very low effective permittivity close to zero, while the other substrate is a standard material with a typical permittivity value.Utilizing an epsilon-near-zero (ENZ) substrate transforms a passive metasurface into an active one.Our approach successfully meets both dual reactive impedances (inductive and capacitive) and impedance featuring a uniform reactive component alongside resistance with positive and negative values (gain and loss).We utilize non-Hermitian impedance boundaries, satisfying parity-time symmetry, to create a planar waveguide that supports LWs.

Results and discussion
Figure 1a shows the schematic of interfaced impedance surfaces that can be considered as a layout designed to support LWs.The presence of two surface modes with perpendicular polarizations can be explained by homogeneous surface impedances in the following manner 6 where η 0 denotes the intrinsic impedance of free space, ζ is a real value, and R, representing the real component of surface impedance, can be either negative or positive.When ζ TM = ζ TE , inductive and capacitive impedance surfaces facilitate the propagation of symmetrical LWs.Since the LW with spin-up can exclusively move forward and the LW with spin-down can only go backward, the spin-momentum locking characteristic enforces the unidirectional nature of these states 26,27 .Indeed, according to the theorem of pseudospin states, two impedance surfaces preserving the EM duality satisfy mirror reflection symmetries.This leads to forming decoupled LWs known as pseudospin states 28 .In an ideal scenario, electromagnetic waves exhibit a singularity precisely at the interface line, as depicted in Fig. 1b.Additionally, surface impedances associated with non-matching electromagnetic responses support a quasi-line mode with a distribution resembling that of the edge modes.
Considering that graphene has the capability to support both TM and TE surface waves, one may utilize freestanding graphene to facilitate the propagation of LWs [29][30][31] .The equations provided below determine the dispersion characteristics of surface waves that find support on a graphene sheet 32 (1) where ω is the angular frequency, c is the speed of light in vacuum and σ g is the surface conductivity of graphene that can be derived using the Kubo formula as 33 Here, we have E f representing the Fermi energy, as the reduced Planck constant, ω as the angular frequency, e as the charge of an electron, T = 300 K for temperature, and τ = 10 −13 s as the relaxation time.Figure 2a shows the real part of graphene conductivity plotted versus frequency at different E f values.The actual value of σ g remains positive across all frequencies and E f values, as depicted in Fig. 2a.Following Eqs. (3),and (4), when σ g is real, and the expression under the radicals is negative, both β TM and β TE become purely imaginary, indicating the propagation of evanescent surface modes.In essence, real values of σ g influence the dissipation loss of waves that are supported by the graphene sheet.
The imaginary component of σ g can change its sign when E f is switched from negative to positive or vice versa, as depicted in Fig. 2b.Consequently, the graphene structure has the capability to support either TE or TM surface waves depending on the value of E f .To illustrate, we consider the propagation of surface plas- mons in graphene at a frequency of 50THz .When E f = 0.4 , the graphene surface conductivity becomes σ g = 6.45 × 10 −6 − j1.39 × 10 −4 , indicating the propagation of TM surface waves in graphene.In this scenario, the normalized propagation constant is Re[β/k 0 ] = 37.94 , indicating strong confinement of the wave to the structure and a guided wavelength for TM surface wave propagation significantly shorter than the wavelength in free space, i.e., guided-TM ≪ 0 .
At E f = 0.1 , the graphene surface conductivity becomes σ g = 3.41 × 10 −5 + j1.95 × 10 −6 .This corresponds to TE surface wave propagation with a normalized propagation constant of Re[β/k 0 ] = 1 .Therefore, in this case, TE mode propagation in graphene closely resembles TEM mode propagation in free space, where the guided wavelength guided-TE is approximately equal to the free-space wavelength 0 .It is important to note that this equality pertains to phase velocity, while field distributions may differ.For the propagation of a LW in a structure, it is essential for the phase velocities of surface waves with TM and TE polarization to be equal.However, freestanding graphene fails to meet this criterion, rendering it incapable of guiding LWs 7 .

Hermitian line waves
Since implementing LWs using freestanding graphene is not feasible, one approach to achieve these waves is by utilizing graphene-based metasurfaces.In the realm of infrared wavelengths, metallic metasurfaces are unsuitable for LW implementation due to the absence of a well-established classical model for metal conductivity in this range.To address this challenge, we propose employing a graphene patch.We present two structures for guiding LWs: a dual-band waveguide that can support LWs in the terahertz and infrared spectra, and a non-Hermitian platform created by placing a graphene-patch metasurface on an ENZ substrate.We provide a numerical analysis (5) of LW propagation in both Hermitian and non-Hermitian metasurfaces.The implementation of LWs has been limited to the terahertz regime.However, our proposed design opens up the possibility of implementing LWs in the infrared range.Non-Hermitian behavior can be achieved by introducing gain and loss elements into the graphene patch, which is placed on the ENZ substrate.Our approach to implementing non-Hermitian LWs is rooted in the parity-time symmetry.Figure 3a illustrates the schematic of our proposed structure, where graphene patch metasurfaces of different unit cell sizes are employed.The utilization of metasurfaces with different sizes aims to create impedances of both inductive and capacitive types at two distinct frequencies by altering E f .To alter the Fermi energy levels of these unconnected graphene patches simultaneously, we propose the use of an extensive set of bias control pads positioned beneath the framework.This allows for the implementation of reconfigurable circuits so that LWs may be guided through optional pathways 7 .Although we consider two different channels for guiding LWs, this concept can be extended to encompass a more comprehensive, multi-channel scenario, allowing LWs to appear at multiple frequencies.
The suggested metasurface comprises regularly arranged subwavelength graphene patches within its structures.Just like other arrangements, these periodic arrays can be represented using equivalent impedance or equivalent surface conductivities in modeling.A closed-form expression is provided for equivalent impedance corresponding to the array of graphene patches which can be written as 34,35 In this context, with p denoting the period, w indicating the patch width, and ε r set at 3.9 as the substrate's dielectric constant, the electromagnetic behavior of the graphene patch can be seen as a series circuit comprising resistance (R), inductance (L), and capacitance (C).Generally speaking, and by way of illustration, we opt for SiO 2 as the substrate layer.The initial component, Z s1 , represents the series R-L resulting from the combination of the graphene surface impedance and geometric parameters.The second part, Z s2 , corresponds to the capacitive response of the graphene impedance surface, which is dependent on the patch's geometry and the background medium's permittivity.It is important to note that the surface impedance of the graphene patch can be adjusted by altering the chemical potential.Nonetheless, there is a restriction on how much the chemical potential of graphene can be raised.To illustrate, when a gate voltage is applied to the graphene sheet through an electrostatic bias, it can elevate the carrier concentration, consequently affecting the values of E f 36,37 .Alternatively, the use of multilayer graphene sheets can lead to achieving higher Fermi energy.As a result, the total conductivity increases proportionally to the number of layers 7,38 .In this context, we examine a four-layer graphene patch metasurface as a means to attain higher levels of E f .This, in turn, leads to a reduction in the real component of the surface impedance of the graphene patch, subsequently improving the propagation characteristics of the LW. Figure 3b illustrates the variation of the imaginary part of the graphene patch impedance in ( 6) www.nature.com/scientificreports/relation to changes in E f for different unit cell sizes at a frequency of 7THz .When p = 10 µ m and w = 8 µ m, it is possible to configure the impedance of the graphene patch to be inductive, whereas for p = 0.3 µ m and w = 0.15 µ m, the graphene patch's impedance can be set as capacitive.Hence, achieving the duality condition for impedances at this frequency becomes readily feasible.Likewise, in Fig. 3c, we depict the imaginary part of the graphene patch impedance at the frequency f = 19 THz .When p = 0.3 µ m and w = 0.15 µ m, the surface impedance can be either inductive or capacitive, contingent upon the appropriate choice of the E f value.Conversely, for p = 10 µ m and w = 8 µ m, the impedance can be tuned to become inductive.Following a similar approach, electromagnetic duality can be readily established at this frequency, enabling the creation of a dual-band line-wave waveguide by manipulating E f values at different frequencies.Figure 3d displays the electric field distribution of the LW at operational frequencies.At a frequency of f = 7THz , we have chosen E f values as follows: E f 1 = 0.56 , E f 2 = E f 3 = 0.109 .At this frequency, the corresponding impedance values are Z 1 = 30 + j188.5 and Z 2 = Z 3 = 277.6 − j754 .Note that E f (i) corresponds to Z (i) .Similarly, at a frequency of f = 19THz , we set E f values to be E f 1 = 1.26 , E f 2 = 0.3 , and E f 3 = 0.99 .The corresponding impedance values are Z 1 = 13.5 + i188.5 , Z 2 = 97.57+ i188.5 , and Z 3 = 27.2 − i754 .The proposed dual-band waveguide supports the propagation of two distinct frequency bands while maintaining a compact size.Furthermore, the line-wave waveguide can be tailored to support frequency bands that either overlap or do not overlap, depending on the specific requirements of the application.

Line waves applications
LWs possess a range of distinctive characteristics that render them highly appealing for diverse applications.Among the notable advantages of line-wave waveguides is their one-dimensional nature, allowing for the creation of exceptionally compact waveguides for efficient energy confinement and transportation.This quality is particularly attractive in the realm of integrated photonics, where the demand for compact and efficient waveguides is pronounced.LWs also excel in mode confinement, rendering them valuable for applications involving interactions between light and matter, as well as chiral quantum processes.They can be effectively harnessed in the design of photonic devices aimed at manipulating the spin or polarization of light.Moreover, the capability to actively control LWs within graphene-based metasurfaces holds the potential for developing reconfigurable integrated circuits.These circuits could be dynamically adjusted to perform various functions, promising greater flexibility and efficiency in electronic devices 6,7 .
LWs present a promising foundation for signal processing applications, including the realization of logic gates operating in the infrared spectrum.The proposed configuration has the potential to be converted into a versatile waveguide by employing a substantial array of bias control pads positioned beneath the structure.This adaptability offers extensive circuit functionality and the capability to employ it for implementing logic gates.www.nature.com/scientificreports/ to redirect the waveguide input toward the desired output.The selection of the desired output is accomplished through addressing, achieved by regulating E f .In fact, the E f values serve as the selection lines for the 1-to-4 demultiplexer.For instance, by controlling the impedance surface values, the input can be directed to output 3. Similarly, this configuration can be utilized to create a NOT gate.Since purely ohmic impedance surfaces lack the capacity to support surface waves, they can serve as a filtering component.This sets the stage for the implementation of a NOT gate, as depicted in Fig. 4b.This functionality hinges on the alteration of impedance surfaces through the adjustment of E f .By modifying E f , it becomes feasible to obstruct the wave's propagation towards the output path.Consequently, the structure functions as an inverter, converting a logic level of 1 to 0. Conversely, to convert a logic level of 0 to 1, a control signal can be employed.In this scenario, when the waveguide input is 0, the output can be set to the logic level of 1 39 .
It is worth noting that, similar to gradient metasurfaces, graphene patches offer the potential to finely manipulate the phase, intensity, and propagation characteristics of electromagnetic fields.Specifically, this enables precise phase modulation of electric fields in response to ζ variations, as depicted in Fig. 4c.This capability makes it feasible to design a straightforward phase shifter without the need for added complexity. Figure 4d illustrates how the electric field distribution of the LW varies at different ζ values.The LW exhibits greater field intensity at lower ζ values.Furthermore, as ζ increases, the guided wavelength of LWs also increases, whereas LWs with lower ζ values remain tightly confined to the interface line.This adjustability in the LW characteristics adds to its appeal for radiation-related applications.It should be noted that the distribution of the fields is calculated at the frequency of 7 THz.

Non-Hermitian line waves
LWs have been observed when two complementary impedance surfaces converge to form a one-dimensional channel.In this arrangement, one of the surfaces must exhibit inductive characteristics, facilitating the propagation of transverse magnetic surface waves.Conversely, the other surface should display capacitive properties, enabling the support of transverse electric surface waves.Moreover, research has indicated that if non-Hermitian metasurfaces adhere to parity-time symmetry conditions, line waves emerge at the interface of non-complementary surfaces.Consequently, LWs can materialize at the junction of impedance surfaces that are exclusively inductive or capacitive, influenced by a gain and loss mechanism.
It is worth noting that various arrangements of non-Hermitian metasurfaces can be explored to create novel waveguides capable of supporting LWs. Figure 5a illustrates the geometric layout of a junction formed by metasurfaces adhering to parity-time symmetry conditions, enabling the implementation of non-Hermitian LWs.In this setup, both non-Hermitian inductive and capacitive metasurfaces are employed to construct a parallel plate waveguide.The impedances composing the waveguide conform to the following parity-time symmetry condition The impedance surfaces mentioned earlier, which create the waveguide boundaries, may not always exhibit mirror reflection symmetries.As a result, this waveguide does not facilitate the propagation of pseudospin states.Figure 5b shows the electric field distribution of LWs at the frequency of 10 THz.The impedance surfaces are spaced apart by a subwavelength distance, leading to the support of coupled LWs.www.nature.com/scientificreports/

Possible implementation of the non-Hermitian waveguide
A prevalent method for inducing non-Hermitian behavior involves introducing gain and loss effects through strategic positioning of gain materials, lossy dielectrics, or active elements within the meta-atoms [40][41][42][43][44] .These phenomena can be achieved at particular frequencies or across a broad spectrum, contingent upon the metasurface's design.A practical realization of this concept could involve the utilization of photoexcited graphene 45 .Nevertheless, we employ a novel method to attain this outcome.It is important to note that when it comes to graphene-based metasurfaces on conventional substrates, higher chemical potential values are necessary to broaden the LWs' operational range 7 .
To tackle this problem, our solution involves the use of a metasurface made of graphene patches, which demonstrates both capacitive and inductive electromagnetic responses.This metasurface is positioned on a material with epsilon-near-zero (ENZ) properties, as depicted in Fig. 6a.We opt for aluminum-doped zinc oxide as the ENZ substrate, and describe the complex permittivity of this material by applying the Drude oscillator model [46][47][48] where ε b = 3.8 represents the high-frequency permittivity, ω p = 0.5 × 10 15 corresponds to the plasma frequency, which is directly related to the concentration of free carriers, and γ = 0.1 × 10 13 denotes the Drude damping rate.These parameters collectively influence the ENZ frequency at which ε r reaches zero within the infrared range, as illustrated in Fig. 6b.
It is worth noting that the surface impedance of the graphene patch has a unique dual property (that is, it can either be inductive or capacitive).This sets it apart from subwavelength metallic patch/grid metasurfaces, which typically exhibit a predominantly reactive response.This dual characteristic of the graphene patch allows for the manipulation of the electromagnetic response of the structure by adjusting E f while keeping the frequency constant.Figure 6c provides insight into the surface impedance of the graphene patch at a frequency of f = 42.2THz with ε = 0.24 + j0.01 , p = 0 /10 , and w = 4 × 0 /50 in which 0 is the wavelength which corresponds to the operating frequency.When coupled with an ENZ substrate, the EM response of the graphene patch shifts from being capacitive at lower E f values to becoming inductive at higher E f values.( 9) www.nature.com/scientificreports/ The equations below can be employed to compute the dispersion relationships for the equivalent surface impedance of the proposed structure [7][8][9][10][11][12][13][14][15]49 In Fig. 6d, we can observe the confinement factor of surface waves supported by the graphene patch at room temperature. Fohe mentioned parameters, within the range of 0.028 ≤ E f ≤ 0.075 , a TE mode confinement factor ranging from 1 ≤ Re[β/k 0 ] ≤ 39.85 is obtained.However, as E f increases, a TM mode emerges.In this scenario, within the range of 0.096 ≤ E f ≤ 0.84 , a confinement factor of 1.002 ≤ Re[β/k 0 ] ≤ 31.2 is achieved, and it is worth noting that a three-layer graphene patch is utilized on the TM side to compute the dispersion relation of the TM surface wave.The inclusion of a multi-layer graphene patch serves the purpose of minimizing dissipation losses.Additionally, the use of an ENZ substrate enables the implementation of LWs with low values of E f .
In practice, achieving the condition of parity-time symmetry, which involves reversing the resistance, is a considerably difficult challenge.Instead, due to the fact that electrons and holes in optically pumped graphene have energy spectra without energy gaps, negative resistance occurs specifically in the terahertz frequency range 50,51 .Our suggested design opens the door to introducing the gain-loss phenomenon in the infrared spectrum.This facilitates the utilization of LWs by harnessing both duality and parity-time symmetries.Since the material that is used as an ENZ substrate has a complex permittivity, according to Eq. ( 6) the real part of surface impedance becomes negative at some specific frequencies.This means that the structure can be used in an attempt to implement waveguiding systems that support non-Hermitian LWs.Consequently, the proposed structure is a feasible way to attain the non-Hermitian behavior.If a complex-valued permittivity is considered as ǫ r = Re(ǫ r ) + Im(ǫ r ) , the real part of graphene patch surface impedance according to Eq. ( 6) can be written as follows According to Eq. ( 12), Re(Z s ) can be negative providing its first part is smaller than the second part.This can be achieved at some frequencies with a lossy ENZ substrate.Figure 7a  .03for a one-layer metasurface.These values are notably higher compared to the scenario where a conventional substrate is used, and their upper limit is substantially increased.Consequently, this leads to an expanded operational range for LWs by controlling the ζ values.It is critical to bear in mind that the proposed layout can be used to implement both Hermitian and non-Hermitian LWs.(10)   Figure 8b, depicts the confinement factor of the line mode at a frequency of 19 THz.For ζ values ranging from 1.3 to 2, a higher confinement factor within the range of 2.62 ≤ Re[β] ≤ 8.56 is achieved.However, there are limitations on the ζ values.It is worth noting that all limitations pertaining to the dispersion characteristics of the line mode are dictated by the real and imaginary parts of surface impedances.As discussed earlier, noncomplementary metasurfaces support a quasi-line in which a portion of electromagnetic energy is restricted within the supporting metasurfaces.Consequently, this mode exhibits a striking similarity to well-known edge modes as shown in Fig. 8c.Lastly, Fig. 8d

Discussion
To implement non-Hermitian line-wave waveguides in practice, it is crucial that the parity-time symmetry should be satisfied.According to Eqs. ( 7) and ( 8), the parity-time symmetry condition necessitates complete equality in the imaginary part of surface impedances, and the absolute value of the real part of surface impedances should be identical, with one surface having negative resistance and the other positive resistance.As previously discussed, in the case of LWs guided by Hermitian surfaces, inequality in the values of ζ leads to the appearance of the quasi- line wave.Conversely, unequal real parts of impedances in Hermitian LWs do not impact the mode symmetry and do not hinder the propagation of LWs.The practical implementation of Hermitian LWs over a wide range of ζ values can be easily achieved.However, implementing quasi-line waves using non-Hermitian surfaces is not feasible due to the violation of the parity-time symmetry condition.Particularly, this feature holds true for the non-Hermitian inductive impedance surfaces.It is worth noting that unequal resistance values also violate the parity-time symmetry condition.Indeed, precise compliance with this condition results in the appearance of non-Hermitian LWs, imposing significant limitations on the practical implementation of non-Hermitian line modes within specific frequency ranges.For better comparison, Fig. 9a illustrates the distribution of Hermitian LWs guided by lossy dual surfaces, enabling symmetrical line mode guiding.The field distribution, more commonly For a clearer comparison, Fig. 9c showcases the distribution of non-Hermitian LW modes for Z TM1 = 0.5 × η 0 + j0.8 × η 0 , and Z TM2 = −0.5 × η 0 + j0.8 × η 0 , where the parity-time symmetry relation is basically satisfied, resulting in a symmetric field distribution.As previously discussed, the non-adherence to parity-time symmetry hampers the propagation of non-Hermitian LWs.For instance, Fig. 9d displays the propagation of LWs in a non-Hermitian waveguide, where a disturbance is created towards the end of the waveguide by selectively choosing impedance values as Z TM1 = 0.5 × η 0 + j0.8 × η 0 , Z TM2 = −0.5 × η 0 + j0.8 × η 0 , Z TM2 = 0.5 × η 0 + j0.8 × η 0 , and Z TM4 = −0.5 × η 0 + j0.6 × η 0 .Obviously, this disturbance prevents wave propagation.It is important to note that applying this disturbance to the real part of impedances would also inhibit the appearance of non-Hermitian LWs.Therefore, practical limitations arise concerning non-Hermitian LWs when considering the complex requirement of parity-time symmetry.
Our proposed waveguiding structure offers several distinctive advantages over previously studied structures.For instance, the proposed Hermitian line-wave waveguide supports one-way modes called "pseudospin states" 28 .Table 1 presents a comprehensive comparison between the characteristics of the proposed waveguides and the structural configurations investigated in the literature (references 7,14 .This table meticulously outlines and contrasts various structural features, offering a detailed overview and analysis of the similarities and differences observed among these different waveguide designs.In the tabe, "N/A" stands for "Not Applicable".

Conclusion
In summary, we introduce a graphene patch with a unique dual capacitive-inductive property, offering the potential for both Hermitian and non-Hermitian line-wave waveguides.We explore various applications for these waveguides and demonstrate that placing the graphene patch on an ENZ substrate facilitates the creation of waveguides with low E f values.By adjusting E f and the number of layers in the multilayer graphene structure, we can control both positive and negative resistance values in the structure, thereby influencing the transition characteristics of line waves.We also highlight how modifications to ζ can impact line wave intensity, phase, and confinement, suggesting possibilities for reconfigurable circuit development and controllable logic gates.Furthermore, we introduce different planar waveguide configurations capable of supporting line waves, including a dual-band waveguide, and a parallel plate waveguide that guides coupled line waves.

0 2 2 Figure 1 .
Figure 1.Unique characteristics of the LW.(a) A one-dimensional junction formed by interfacing dual impedance boundaries, and (b) singularity of the LW along the interface line of dual impedance surfaces.The LW decays exponentially in both transverse directions.

Figure 2 .
Figure 2. Graphene surface conductivity.(a) Real, and (b) imaginary parts of the graphene conductivity versus frequency and E f variations at room temperature.

Figure 3 .
Figure 3.The proposed structure and its related characteristics.(a) A schematic that illustrates the proposed dual-band line-wave waveguide that utilizes a graphene-based metasurface, (b) the graph depicting the surface reactance of the structure at the frequency of ω/2π = 7THz versus chemical potential, (c) surface impedance of the graphene patch for different unit cell sizes at the frequency of ω/2π = 19THz , and (d) electric field magnitude distribution at two different frequencies on top of the waveguide.

Figure 4 .
Figure 4. Illustration of possible applications of the LW.(a) A possible implementation of optical circuits such as a demultiplexer, and (b) a Not gate, (c) phase shift of the LW that corresponds to surface impedance varations across the interface line, and (d) control of the intensity and guided wavelength of the LW using ζ variation.

( 7 )Figure 5 .
Figure 5. Non-Hermitian LWs.(a) The geometry of the problem, which involves a non-Hermitian impedance junction that can support coupled line waves, and (b) the distribution of the magnitude of the electric field of the LW at the cross-section of the waveguide.

Figure 6 .
Figure 6.Characteristics and details of the proposed dual impedance surface.(a) Schematic of the proposed dual metasurface on an ENZ substrate and associated circuit model, (b) real and imaginary parts of aluminiumdoped zinc oxide permittivity as a function of frequency, (c) the imaginary part of graphene patch impedance as a function of chemical potential at f = 42.2THz , and (d) the associated normalized propagation constant of surface waves supported by the proposed structure.
illustrates the resistive component of the proposed impedance surface.A shift from positive resistance (indicating loss) to negative resistance (indicating gain) occurs at ε = 0.05 +j0.014 and a frequency of 41THz .In Fig. 7b, one can see the ζ values plotted against E f .The maximum values of ζ are achieved with ζ TM ≤ 8.43 for a three-layer graphene metasurface and ζ TE ≤ 7

Figure 7 .
Figure 7. EM response of the graphene patch positioned on the ENZ substrate.(a) Gain effect achieved by properly adjusting of the E f value, and (b) ζ values of surface waves versus E f .
Figure8aillustrates the dispersion characteristics of the line wave in the dual-band waveguide at a frequency of 7 THz.To analyze the dispersion features of the line mode, we consider impedance surfaces with complex values.It is important to note that the ζ range at each frequency is determined by the values of the surface impedances.For ζ values ranging from 1.6 to 6.48, a confinement factor within the range of 1.18 ≤ Re[β] ≤ 3.42 is achieved.The confinement factor decreases proportionally with an increase in ζ , while for ζ ≤ 1.3 , there is an evanescent wave along the propagation direction, indicating that Re[β] ≤ Im[β] .Conversely, ζ < 1 corresponds to the LWs cutoff, in accordance with observations.Figure8b, depicts the confinement factor of the line mode at a frequency of 19 THz.For ζ values ranging from 1.3 to 2, a higher confinement factor within the range of 2.62 ≤ Re[β] ≤ 8.56 is achieved.However, there are limitations on the ζ values.It is worth noting that all limitations pertaining to the dispersion characteristics of the line mode are dictated by the real and imaginary parts of surface impedances.As discussed earlier, noncomplementary metasurfaces support a quasi-line in which a portion of electromagnetic energy is restricted within the supporting metasurfaces.Consequently, this mode exhibits a striking similarity to well-known edge modes as shown in Fig.8c.Lastly, Fig.8ddemonstrates the propagation characteristics of a quasi-line mode.This mode exhibits a confinement factor ranging from 1.05 ≤ Re[β]/k 0 ≤ 9.4 for ζ values ranging from 0.2 to 50.Unlike the conventional line mode, a quasi-line mode exists for ζ < 1 , provided that ζ TM × ζ TE > 1.

Figure 8 .
Figure 8. Eigen-mode analysis.(a) Dispersion characteristics of the line mode for the dual-band waveguide at the frequency of f= 7THz, (b) at the frequency of f = 19 THz, (c) E-field distribution of the quasi-line mode at the waveguide cross section, and (d) the dispersion diagram of the quasi-line mode.

Figure 9 .
Figure 9. LWs at the interface of Hermitian, and non-Hermitian impedance boundaries.(a) LWs supported by dual impedance surfaces that are lossy.(b) Quasi-line waves at the interface of dual impedance surfaces having unequal ζ values.(c) LWs supported by non-Hermitian inductive impedance surfaces.(d) LWs filtering caused by introducing a defect in the system.

Table 1 .
Comparison of key parameters of reported waveguides.